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ABSTRACT 

In this study, we investigate the shapes of starless and protostellar cores using hydrodynamic, 
self-gravitating adaptive mesh refinement simulations of turbulent molecular clouds. We simulate 
observations of these cores in dust emission, including realistic noise and telescope resolution, and 
compare to the observed core shapes measured in Orion by Nutter & Ward-Thompson (2007). The 
simulations and the observations have generally high statistical similarity, with particularly good 
agreement between simulations and Orion B. Although protostellar cores tend to have semi- major axis 
to semi-minor axis ratios closer to one, the distribution of axis ratios for starless and protostellar cores 
are not significantly different for either the actual observations of Orion or the simulated observations. 
Because of the high level of agreement between the non-magnetic hydrodynamic simulations and 
observation, contrary to a number of previous authors, one cannot infer the presence of magnetic 
fields from core shape distributions. 

Subject headings: ISM: clouds - kinematics and dynamics- stars:formation - methods: numerical — 
hydrodynamics - turbulence 



1. INTRODUCTION 

A successful theory of star formation must explain 
certain basic characteristics of the early stages of cores 
and stellar natal conditions. One property of interest is 
the shape distribution of starless and protostellar cores, 
which is likely related to the the initial conditions of star 
formation such as the local magnetic field configuration, 
level of turbulence, and core collapse timescale (see re- 
view by McKee & Ostriker 2007). Strong magnetic fields 
in the early stages of core formation may either support 
gas perpendicular to the field lines yielding a distribution 
of oblate cores (Mouschovias 1976) or compress the cores 
into a prolate geometry (Fiege & Pudritz 2000) . The for- 
mer argument assumes predominantly poloidal magnetic 
fields, whereas the latter work includes a toroidal field 
component resulting in a helical field geometry. Inde- 
pendent of magnetic fields, prolate cores may also arise 
as an artifact of filamentary cloud geometry, in which 
cores fragment at intervals according to the characteris- 
tic Jeans length of a cylinder (Hartmann 2002). 

Recent large scale turbulent simulations with and with- 
out magnetic fields tend to find distributions of cores that 
contain predominantly triaxial cores when viewed in 3D 
(Klessen k Burkert 2000; Gammie et al. 2003; Basu & 
Ciolek 2004; Li et al. 2004; Offner et al. 2008a). In con- 
trast, simplified numerical studies including ambipolar 
diffusion tend to find triaxial cores with an inclination 
towards oblateness (Basu & Ciolek 2004; Ciolek & Basu 
2006). Observations of molecular cloud cores in various 
star forming regions tend to find projected core aspect 
ratios, q — a/h^ around 2:1 (Myers et al. 1991; Jijina et al. 
1999; Nutter & Ward-Thompson 2007), but comparison 
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of these observations with simulations is hampered by the 
projection of the observed cores onto the plane of the sky. 
Some authors have presented analytic work attempting 
to overcome this difficulty by "de-projecting" the ob- 
served cores statistically. For example, Ryden (1996), 
assuming axisymmetry, finds that cores are significantly 
more likely to be randomly oriented prolate objects than 
oblate objects. Tassis (2007, henceforth T07) utilizes a 
maximum-likelihood method to generate a distribution of 
ellipsoid axial ratios. He uses two base probability dis- 
tribution functions and finds that oblate or triaxial cores 
agree well with observations of cores in Orion (Nutter 
& Ward-Thompson 2007, hereafter NWT), results which 
are insensitive to the assumed underlying distribution. 
Overall, T07 finds that prolate cores are rare, and his 
method rules out a uniform distribution of oblate, pro- 
late, and triaxial cores with greater than 99% confidence. 

Our study is complementary to this previous work, 
but we perform the core shape comparison in the ob- 
servational domain rather than the theoretical one. This 
has a significant advantage over alternative approaches, 
because it allows us to realistically simulate the effects 
of finite telescope resolution and sensitivity, and to re- 
duce the simulated observations and fit core shapes using 
the same methods used for the real data. In this paper, 
we post-process simulations in this manner and compare 
with the dataset of Orion reported by NWT, who collate 
and reanalyze various SCUBA observations of the Orion 
A and B North and South molecular cloud complexes. 
They report the masses and sizes for 393 cores. Con- 
sequently, this dataset not only concerns an interesting 
and well-studied region, but it is also sufficiently large 
for meaningful statistics. 

In section 2, we describe the details of the simulations 
and post-processing. We consider simulations where tur- 
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bulence is continually driven and where it is allowed to 
decay. In section 3, wc present comparisons between the 
simulations and observations and compare with T07. In 
section 4, we summarize our conclusions. 

2. SIMULATED OBSERVATIONS 

As described in OfFner et al. (2008a, b). our two siirm- 
lations are periodic boxes containing an isothermal, non- 
magnetized gas that is initially not self-gravitating. After 
driving turbulent motions in the gas for two box cross- 
ing times, self-gravity is turned on. In one simulation 
energy injection is halted and the turbulence gradually 
decays, while in the other turbulent driving is maintained 
so that the cloud satisfies energy equipartition. A sink 
particle is introduced when the Jeans conditions is ex- 
ceeded on the finest AMR level, where the cell spacing 
is Aa:;=200 AU (Krumholz ct al. 2004). Since isothermal 
self- gravitating gas is scale free, it is easy to normalize the 
simulations to the conditions observed in Orion using the 
thermal Jeans length and thermal Jeans mass (see scal- 
ing relations in Offner et al. 2008b). For the simulation 
normalization, we adopt a gas temperature of T = 20 K 
(NWT). We choose a gas density of p = 9.74 x lO^^i g 
cm^ (number density of hydrogen nuclei uh = 4.2 x 10"^ 
cm~^), corresponding to a simulation box length of 2 
pc. At this density, the typical box column density is 
Nh = 2.6 X 10^^ cm^^. close to the measured central col- 
umn density of the Orion B cloud, Nh = 2.8 x 10^^ cm~^ 
(Maddalenaet al. 1986; Johnstone et al. 2001). Although 
Maddalena et al. (1986) observe a much larger region 
than the individual north and south complexes in Orion 
A and B, column densities are observed to be roughly 
scale-independent (Larson 1981; Heyer et al. 2008) (al- 
though it has been argued that this is a selection effect - 
e.g. Ballesteros-Paredes 2002). Thus, our normalization 
is consistent. 

We pause at this point to add a caveat about our nor- 
malization; our box size is smaller than the individual 
Orion A and B North and South complexes, which have 
projected sizes of ~ 3 — 4 pc. Equivalently, our simu- 
lation contains a smaller number of total thermal Jeans 
masses of gas than the entire Orion A and B North and 
South complexes. The justification for this is that the 
periodic geometry of the simulation enables us to model 
a piece of the cloud rather than the entire cloud. As 
long as we have the correct mean density, the behavior 
of structures that arc much smaller than the box size 
such as cores should not depend on the size of the sim- 
ulation box. Indeed, the similarity between the results 
for driven and decaying turbulence, and the insensitivity 
of our results to our assumed density normalization (see 
below), seem to support the hypothesis that small-scale 
structure is insensitive to such large-scale features as the 
total box size and the number of thermal Jeans masses 
it contains. However, establishing this point definitively 
would require a suite of simulations with varying total 
sizes but the same resolution, and such a study is un- 
fortunately too computationally costly to perform. We 
therefore simply caution readers on this point and pro- 
ceed. 

We run the simulations with gravity for a global freefall 
time, tfi, and we compare with the observations at f4[/2 
and tg. The column density of the two simulations at tg 
is displayed in Figure 1. At these times, the Mach num- 



ber, A4, is 8.4 for the driven simulation and Ai ~ 5.3, 4.5 
for the decaying simulation. For the driven simulation, 
this corresponds to half the Mach number of the larger 
Orion A and B regions that we compare with (Maddalena 
et al. 1986). However, it is unlikely that the high Mach 
number flow regions will have a significant effect on the 
details of the cores, which are generally subsonic to tran- 
sonic. 

It is possible to convert the simulation column density 
to an observed intensity using the relation 

I,^ = Ai'Hf^mbMH'«7/B,y(rcloud) (1) 

where 1^ is the flux density per beam at frequency i^, ^Imb 
is the solid angle subtended by the beam, hh is the mean 
mass per H atom, Ki, is the dust opacity at frequency v, 
and B^(T) is the Planck function (Enoch et al. 2007). 
We set £^=850 fim to match the observations of NWT, 
and following them we adopt ksso = 0.01 cm^g~^. Note 
that eqn. 1 assumes that the gas is optically thin at 850 
fim. The densest sightlinc through the simulation has a 
column density of E = 0.6 g cm~^ or T85o=0.06, which 
is safely optically thin. 

In order to facilitate comparison with NWT, we post- 
process the simulation data to have noise and resolution 
comparable to the SCUBA data set. In each projec- 
tion direction, we integrate along the line of sight and 
convolve the column density image with a beam of res- 
olution 14" . We assume that the simulated cloud lies 
at a distance of 400 pc, which is the average distance 
adopted by NWT. To each pixel in the smoothed im- 
age, we add a Gaussian noise distribution with ctnwt = 
20 mJy beam~^, correlated over the FWHM size of the 
beam. This reproduces the coarser pixel resolution noise 
inherent in the SCUBA data. 

Having generated a simulated column density map, we 
analyze the data using the same procedure outlined in 
NWT, which we describe in the following steps. First, 
from the post-processed simulation data, wc generate a 
large scale structure map by convolving the data with 
a beam size of 1'. Second, we subtract this map from 
the high resolution version to remove large scale struc- 
ture to make core identification easier. Third, we define 
cores as density enhancements within 3o'nwt contours 
that contain a peak above 5crNWT- With few exceptions, 
core shapes are generally elliptical. Finally, we fit an 
ellipse around each peak and match to the Sctnwt con- 
tours "by eye" guided by a sample eye-fit from NWT 
(D. Nutter, 2008, priv. comm.). An automated fitting 
algorithm would clearly be preferable to the "by eye" 
procedure from the standpoint of reproducible compar- 
ison between different samples and authors. However, 
after experimenting with several options we were unable 
to find an algorithm that adequately reproduced the fits 
from NWT, while we were able to do so reasonably well 
by eye. Since the goal of this paper is to compare to the 
observed sample, we use by eye fits. Based on the ellip- 
tical fits we assign each core an axis ratio q. We identify 
cores that contain a sink particle within 0.05 pc of their 
core center as "protostcllar cores" , while cores that do 
not contain a sink particle are "starless" . In the sample 
of Orion cores, the protostellar cores are identified with 
the Spitzer IRAC camera. 

Figure 2 shows one projection through the simulation 
domain for the driven run, along with the bounding el- 
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lipses for all the cores we identify in that projection. For 
comparison, the unprocessed driven column density is 
shown in the right panel of Figure 1. To improve the 
statistics, we include cores from maps in each cardinal di- 
rection of the simulation. The Kolmogov-Smirnov (KS) 
test provides a good statistical measure of the agreement 
of two distributions (Press 1992). It derives 1 minus the 
confidence level at which the null hypothesis that the two 
were drawn from the same underlying distribution can be 
ruled out^. . In the remainder of the paper, we report 
the quantity 1 minus the confidence level multiplied by 
a factor of a 100 to give a percentage. Using a KS test, 
we find that core samples from different lines of sight 
are consistent with being drawn from the same distribu- 
tion (KS statistics ~ 50%). This implies that the core 
ratios in the three projections are statistically indistin- 
guishable, so our procedure of treating each orientation 
as an independent sample is consistent. The agreement 
among the orientations also shows that fitting by eye is 
reasonably reproducible. 

We verify that the distribution of core ellipses is not 
strongly dependent on the details of the normalization 
by comparing the ellipses in the fiducial case with core 
samples assuming a = O.Sunwt, and a = 2crNWT- This 
is equivalent to adopting the same contour level while 
changing the average density of the simulation by a fac- 
tor of 2. Figure 3 shows the fiducial distribution of axes 
ratios and the two distributions with different contours. 
Wc find that the new axis ratio distributions are consis- 
tent with being drawn from the same distribution at 99% 
and 17% confidence, respectively, when comparing sam- 
ples of equal number. Using a lower or higher contour 
level respectively increase or decreases, respectively, the 
number of cores in the population by ~30%. 

3. DATA COMPARISON 

In table 1 we report a summary of the statistical prop- 
erties of our cores, including the size of the sample, the 
absolute sizes of cores, and their aspect ratios. We com- 
pute these properties for all cores, for the real Orion data 
set, and for both the driven and undriven simulations at 
our two sample times t = ^ts and t = ts. In each case 
we report properties both for the entire population of 
cores regardless of whether they contain stars and for 
the starless and protostellar populations separately. 

3.1. Core Sizes 

^ Formally, the two-sided K-S statistic we use is computed as 
follows. Consider two sets of N and M measurements of some 
quantity x (in our case x is the axis ratio, and the two sets 
of measurements are the simulated and observed values). Let 

Fi^('x) and Gm(^) be the cumulative distribution functions for 
those measurements, i.e. Fj^{x) is the fraction of the N measure- 
ments that yield a value less than or equal to x, and similarly 
for Gm{^)- The Kolmogorov-Smirnov statistic is then defined as 
Dn,m = supj,|_Fjv(^) — Gm(^)\, i-c. the maximum distance be- 
tween the two cumulative distribution functions. Kolmogorov's 
Theorem then states that we can reject the null hypothesis that 
F^(x) and Gm(^) were drawn from the same parent distribution 
with a confidence level a if ^NM/(N + M)Dn.m > Ka, where 
Ka is defined implicitly by the equation 1 — a = Pt{K < Ka) 

and Pr{K < x) = 1 - 2 (^1)'"^'^"^'^''^ • Intuitively a gives 
the probability that we could have measured a value of Djq^M as 
large as we did if Fi^{x) and Gm{^) were actually drawn from the 
same parent distribution. The smaller a is, the less likely that our 
samples would have produced such a large Dm^m if the samples 
were drawn from the same parent distribution. 



In this section, we compare the physical sizes of the 
simulated and observed cores, where the lengths of the 
semi-major and semi-minor axes are given in Table 1. 
Generally, we find that the medians of the net distribu- 
tion of simulated core sizes are ~ 20% smaller. Whereas 
the simulated starless core sizes are fairly similar to ob- 
servation, the observed protostellar core sizes are much 
larger than their simulated counterparts. This discrep- 
ancy is likely an artifact of the sink particle accretion al- 
gorithm: once a core forms a sink particle the surround- 
ing gas is accreted more quickly and to higher masses 
without losses from outflows, so that the reservoir of 
bound gas around the sink particle is rapidly depleted. 
It is also possible that magnetic fields, which we neglect, 
play a significant role in supporting the protostellar en- 
velopes, thus slowing the collapse process and contribut- 
ing to the larger sizes of the observed protostellar cores. 

As shown in figure 4, the simulation core size distri- 
butions have smaller dispersions than the observation. 
Although some of the difference can be attributed to the 
smaller sizes of protostellar cores, the smaller core sam- 
ple and simulation domain size may also contribute. The 
minimum core size is most likely set by the observation 
resolution. 

We can compensate for the rapid depletion of gas in the 
outer parts of protostellar cores if we adopt a significantly 
lower estimate of the telescope noise (cr ~ 0.2c7nwt) and 
thus a lower contour threshold for defining cores. This 
has the effect of making the protostellar cores somewhat 
larger while leaving the starless sizes mostly unchanged, 
so that the overall size distribution is in better agreement 
with the observations. Using this lower noise level does 
not significantly alter the core shape distribution, how- 
ever, which suggests that the discrepancy in core sizes 
is not significant for the purpose of determining core 
shapes. We therefore proceed with our analysis using 
the real telescope noise level, a = ctnwt- 

3.2. Overall Shape Distributions 

As shown in Table 1, we find similar means and medi- 
ans for the shape distributions. The characteristic mean 
falls around q = 0.6 — 0.7. The maximum aspect ratio 
is also very similar in all the cases, around 4:1, and the 
most elongated core is starless. In the simulations, core 
elongation is a result of the initial filamentary gas struc- 
ture out of which the cores form. It is the remnant of the 
turbulence rather than a signature of magnetic fields. 

We next characterize the similarity of the distribu- 
tions by using a KS test. One interesting aspect of 
the Orion dataset is that we can rule out the possibil- 
ity that the separate Orion A and B populations origi- 
nate from the same parent population with > 95% con- 
fidence. This disagreement is caused by the higher frac- 
tion of elongated cores in Orion B. The physical differ- 
ence between the North and South regions are otherwise 
not large: both are sites of high-mass star formation, 
both have similar average column densities (Maddalena 
et al. 1986), and the patches surveyed are roughly the 
same size. However, Orion A has a 25% larger velocity 
dispersion, 5.1 km s~^ as measured in CO (Maddalena 
et al. 1986). The magnitude of this difference is a useful 
number to bear in mind when characterizing the extent 
of agreement or disagreement with the simulations. In 
comparison, the driven and decaying samples are consis- 
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Fig. 1. — Logarithm of the column density, S, for one projection of the decaying (left) and driven(right) simulations at Itff. 
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Core axis ratio b/a minimum, median, and mean and median core sizes. 
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^Observed Orion molecular cloud cores (NWT) 
''Driven turbulence simulation at Itg 

Driven turbulence simulation at ^tg 

Undriven turbulence simulation at Itff 
''Undriven turbulence simulation at i% 

'Median projected semi-major (a) and semi-minor (b) size in units of 100 AU. 



tent with being drawn from the same parent population 
at 88% and 66% confidence for \ts and |iff , respectively. 

We note that the KS statistic is somewhat influenced 
by the size of the distribution being compared. For exam- 
ple the smaller the samples, the more likely the test will 
conclude two samples are consistent with being drawn 
from the same distribution. In general we find if we 
always compare distributions of the same size (by ran- 
domly selecting cores from the larger distribution), KS 
agreement rises by 5 — 10%. 

Table 2 shows the KS statistics for comparisons of the 
axis ratios in Orion to each of the simulations. Overall, 
the decaying turbulence simulation agrees better with the 



Orion populations, although both simulations have fairly 
high agreement with Orion B. Given the uncertainties 
in the normalization and the disagreement between the 
Orion A and B samples, we consider agreement greater 
than 10% to be encouraging. Interestingly, the disagree- 
ment of the driven population is determined mainly by 
the disagreement of ellipses in one particular projection 
at Itff, which has an overabundance of elongated cores. 
Minus the cores in this projection, the driven sample 
agrees with ~ 17% confidence at this time. Figure 5 
shows the cumulative distribution function of the core 
shapes for the observations and two simulations. 
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Fig. 2. — Column / flux density for one projection of the post- 
processed driven simulated cores at 1%. The image includes noise 
and beam-smearing. The 3(tnx contours are marked in black and 
the fitted ellipses are overlaid in white. 



Fig. 3. — Number of cores as a function of axis ratio q = b/a 
along a single projection for cores defined using the fiducial crisiT, 
ScNTi and O.Sctnt at Itff. 




50 100 150 50 100 150 50 100 150 
a [100 AU] a [100 AU] □ [100 AU] 

Fig. 4. — The sizes of major axis, a, in units of 100 AU for the 
total, starless, and protostellar cores, from left to right, at Ug. 



3.3. Starless and Protostellar Core Shapes 

KS agreement for starless and protostellar cores is com- 
pared in Table 3. We find similar agreement with the 
Orion data for both the driven and undriven starless core 
shapes at lift. In contrast with the observational data, 
which has a high level of agreement between the Orion 
starless and Orion protostellar core shapes 60%) the 
driven starless and driven protostellar core shapes are 
only moderately similar (^13%), while the decaying star- 
less and decaying protostellar core shapes are quite dis- 




co 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.5 0.8 1.0 
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Fig. 5. — Cumulative distribution function of the total, starless, 
and protostellar shape distribution, from left to right, at Itfj. 
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KS AGREEMENT OF SIMULATIONS WITH ORION FOR NET 
POPULATIONS 
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^Observed Orion molecular cloud cores (NWT) 
''Driven turbulence simulation 
'^Undriven turbulence simulation 

similar 0.02%). This is illustrated in second and third 
panels of Figure 5, which show the undriven starless dis- 
tribution of shapes significantly to the left of the Orion 
distribution, while the undriven protostellar distribution 
falls to the right. As a result, the individual simulation 
starless and protostellar shape distributions can be much 
less similar to the observed cores than the net simula- 
tion shape distribution (see Ui in Table 2 and 3). One 
caveat of this comparison is that the actual goodness of 
agreement depends not only upon agreement between the 
net shape distributions but upon the agreement between 
the individual starless and protostellar shapes and ob- 
servations. Figure 5 suggests that, particularly for the 
decaying case, the K-S test may overestimate the over- 
all similarity to observed cores. Observational data bears 
out the similarity of the distributions of starless and pro- 
tostellar shapes (Myers et al. 1991; Jijina et al. 1999). In 
the decaying simulation, the difference is likely due to the 
more rounded protostellar cores, which are experiencing 
strong collapse (Offner et al. 2008b). 

In Figure 6 we plot the number of starless cores as a 
function of axis ratio, q. For comparison, we also plot 
the T07 maximum-likelihood curves drawn from normal 
and beta distributions. In spite of the larger variation in 
the simulation data (due to the smaller sample size), the 
simulation data also appears to follow the curves reason- 
ably well. Offner et al. (2008a) report that the proto- 
stellar cores in both simulations are mainly trixial, with 
some preference for prolateness over oblateness. To in- 
vestigate the 3D shape distribution of the starless cores, 
which were not examined in Offner et al. (2008a) , we first 
triangulate the 2D projected positions to identify the 3D 
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Fig. 6. — Histrogram of starless cores ratios q = b/a at Itjf. 
Results from T07 derived assuming underlying beta and normal 
distribution for core ratios have been overlaid. The error bars 
shown on each bin reflect \/N counting statistics. The bins for 
each sample are centered at the same values of q, but the plotted 
points have been offset slightly to the left or right to allow the error 
bars to be distinguished. 



TABLE 3 

KS AGREEMENT OF SIMULATIONS WITH OrION FOR STARLESS AND 
PROTOSTELLAR POPULATIONS 



O'' Starless (%) O Protostellar {%) 



Di'' Starless 1.2 Di Proto 20.7 

Di Starless 2.2x10"^ Di Proto 2.5 

Starless 1.2 Proto 2.9 

Ui Starless 0.04 Ui Proto 0.4 

2 2 



^Observed Orion molecular cloud cores NWT 
''Driven turbulence simulation 
'^Undriven turbulence simulation 

coordinates of the core center. By setting a minimum 
density cutoff for cells within 0.1 pc of the core center, 
we define the gas contained in the core. We adopt a den- 
sity cutoff of the minimum of nn = 2 x 10''' cm~'^ and 
0.2 riHpGak, where protostellar cores generally satisfy the 
former and starless cores the latter. Unlike the core def- 
inition in Offner et al. (2008a), we do not require that 
the gas be bound. We apply principle component analy- 
sis to the set of cells comprising each core to identify the 
eigenvalues of the cardinal axes ( JoUiffe 2002) . As shown 
in Figure 7, the cores for both simulations are predom- 
inately triaxial. The remaining cores are preferentially 
prolate. However, the ratio of the number of prolate to 
oblate cores in 3D is somewhat sensitive to the chosen 
cutoff density. These results appears to be inconsistent 
with the claim of T07 that the observed core axis dis- 
tribution implies that prolate cores are rare while oblate 
cores are more common. 

3.4. Time Dependence of Core Shapes 

We find less KS agreement between the Orion core 
shapes and the simulated cores at ^iff. The origin of 
the disagreement is mainly due to the smaller mean axes 
ratios in both the driven and decaying cases as illustrated 
in Table 1. The difference between starless cores and pro- 
tostellar cores is also more pronounced at this time since 
the starless cores in both simulations are more elongated 
and the agreement between the observed and simulated 
starless distributions is much worse than for the proto- 



FlG. 7. — Plot of the 3D core aspect ratios c/a vs. b/a at Itff. 
The dashed lines indicate the boundary between prolate, triaxial, 
and oblate cores, from left to right. 

stellar core distributions. This may be because the simu- 
lated protostellar cores are necessarily located in regions 
of the flow that are dominated by gravitational effects. 
In contrast, at early times the shapes of starless cores are 
more strongly influenced by turbulence rather than grav- 
ity and so they appear more elongated and filamentary. 
This suggests that gravitational fragmentation, not only 
turbulence, significantly influences core shapes. 

4. CONCLUSIONS 

Overall, we find a high level of similarity between ob- 
servations of core axis ratios in Orion and simulations 
of core formation in a self-gravitating, non-magnetized 
turbulent medium, with either driven or decaying turbu- 
lence with best agreement occurring after one dynamical 
time. The similarity extends from the mean and median 
core axis ratios to the distributions of starless and pro- 
tostellar core shapes. We obtain good agreement despite 
the absence of magnetic fields in the simulations, which 
may indicate that the local magnetic field in Orion is not 
strongly influencing the core shapes. This is supported 
by both turbulent magnetic simulations such as those 
by (Ballesteros-Paredes 2002) and observations that find 
the shapes of high density structures are not strongly 
correlated with the magnetic field direction (Ballesteros- 
Paredes 2007). Moreover, we find that a population of 
cores that is intrinsically triaxial, but with a tendency 
to be more prolate than oblate, can produce an observed 
distribution of core axis ratios consistent with what is 
seen in Orion. 

The axis ratio distributions are also quite similar in 
the simulations with driven and decaying turbulence. In 
fact, we find that the shape distributions of Orion A and 
Orion B are often more dissimilar (to 95% confidence) to 
each other than to the simulations. This indicates that 
the level of turbulence does not play a signiflcant role in 
determining core shapes. Both simulations compare with 
larger confidence to the Orion B core sample, although in 
terms of statistics, the decaying simulation gets slightly 
better agreement with the total Orion sample of cores. 

Increasingly large and complete observational data sets 
invite important comparisons with simulations, which 
could shed light on both the theory of star formation 
and details of the molecular clouds we observe. Our re- 
sults provide a cautionary note that such comparisons 
should preferably be done by projecting from simulations 
into the observational domain, including realistic sensi- 
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tivities and resolutions. Further simulations coneerning 
the effects of magnetic fields, combined with more de- 
tailed simulated observations, would be a beneficial fu- 
ture direction of research. 
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